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Abstract. Frustrated magnets in high magnetic field have a long history of offering beautiful 
surprises to the patient investigator. Here we present the results of extensive classical Monte 
Carlo simulations of a variety of models of two dimensional magnets in magnetic field, together 
with complementary spin wave analysis. Striking results include (i) a massively enhanced 
magnetocaloric effect in antiferromagnets bordering on ferromagnetic order, (ii) a route to an 
m = 1/3 magnetization plateau on a square lattice, and (iii) a cascade of phase transitions in a 
simple model of AgNi02. 



1. Introduction 

It is precisely the feature that makes frustrated magnets so compelhng — the multitude of 
groundstates available — that also makes them so difficult to study, especially if we take in 
account the extra contribution of an external magnetic field. In this article we present a snapshot 
of ongoing work on several experimentally motivated models of two-dimensional frustrated 
magnets in high magnetic field. Our main tool will be classical Monte Carlo simulation, but 
where the standard local update Metropolis algorithm is augmented by a parallel tempering 
scheme [1], which greatly improves the ergodicity of the method. Simulation results will also be 
compared with simple mean field theories and classical spinwave calculations. 

One concept that will be used throughout this work is the magnetocaloric effect (MCE), 
Omce = 1 5) the rate of change of temperature of the system with applied magnetic field 
at constant entropy. This quantity is useful as a tool for identifying phase transitions in both 
simulation and experiment, and forms the basis for magnetic refrigeration. 

2. The J1-J2 Heisenberg model on a square lattice 

The first model we consider is the Heisenberg model on a square lattice with additional 
2"'^-neighbour (J2) bonds 




(1) 



Despite its apparent simplicity, this is a very interesting model, and is thought to describe 
two new classes of oxide compounds [21 [3] . The classical phase diagram in zero field is shown 
in Fig. [H The quantum phase diagram and finite field properties of the spin-1/2 quantum 
model have been extensively investigated for both antiferromagnetic (AF) [1] and ferromagnetic 
(FM) [H El [71 E] Ji. New quantum phases are found in the highly frustrated regions where 
|Ji| = 2J2. Here we explore the role of thermal fluctuations and finite magnetic field for these 
parameter sets. 
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Figure 1. (Colour online) Clas- 
sical phase diagram of the Ji- 
J2 Heisenberg model. Parameters 
for quasi-2D vanadates are taken 
from [7j. 




Figure 2. (Colour online) MCE of the classical 
J1-J2 Heisenberg model for Ji = 1, J2 = 0.5 (T = 
0.01), showing successive phase transitions. T 
and h are measured in units of Ji. The MCE is 
normalised to that of an ideal paramagnet. 



Considering first AF Ji > 0, in Fig.[2]we present the MCE for Ji = 1, J2 = 0.5, as determined 
by classical MC simulation for a cluster of = 24 x 24 spins at T = 0.01. Binned data from 
~ 2 X 10^ MC steps of simulation, preceded by ~ 6 x 10^ MC steps of thermalization, were 
analysed using a jackknife procedure. Analysis reveals four distinct phases, the most striking 
feature being a m = 1/2 plateau (c.f. [9]). 

The MCE also shows interesting structure within the collinear AF phase for FM Ji < 0, 
approaching the classical critical point Ji = — 2 J2, as shown in Fig. [3l Here the high degeneracy 
in the spin wave spectrum shows up as a massive enhancement in the MCE at low fields, also 
discussed for quantum spins in [7] . 
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Figure 3. (Colour online). MCE of the J1-J2 
Heisenberg model on a square lattice for different 
ratios of J2/J1 {Ji < 0), showing a massive 
enhancement near the transition from collinear 
(7r,0) Neel AF to FM at J2/J1 = -0.5.Arrows 
show rough values of the saturation field in each 
case. Results are for 24 x 24 spins, at T = 0.01, 
normalised to the maximum value of the MCE in 
a nearest-neighbour Heisenberg AF. T and h are 
measured in units of Ji. 



3. The J1-J2-J3 Heisenberg model on a square lattice 

A new series of layered Dion-Jacobi compounds containing spin-1/2 Cu{3dP) ions on a square 
lattice have recently been discovered [3]. These materials are believed to have FM nearest 
neighbour exchange Ji, but not all of their properties can be reconciled with a simple J1-J2 
model. In particular (CuBr)Sr2Nb30io exhibits an unexpected m = 1/3 magnetization 
plateau [10]. This has motivated us to consider the quantum and classical properties of Eq. [T] 
with additional AF third neighbour interaction J3. 

So far as classical MC simulations are concerned, our main results are summarised in Fig. HI A 
clear m = 1/3 plateau is found for Ji = —1, J2 = 1.0, J3 = 0.5. The classical h-T phase diagram 
clearly has much in common with that for the nearest-neighbour Heisenberg AF on a triangular 
lattice [TT]. The m = 1/3 plateau is also observed in exact diagonalization calculations for the 
equivalent spin-1/2 J1-J2-J3 model [12j. These results clearly suggest that this model deserves 
further study, both as a "toy model" for (CuBr)Sr2Nb30io and as an interesting problem in its 
own right. 
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Figure 4. (Colour online.) Left Panel : Magnetic phase diagram for the classical J1-J2-J3 
Heisenberg model on a 24 x 24 site square lattice, for Ji = —1, J2 = 1, J3 = 0.5. Points are 
taken from simulation results for magnetic susceptibility x, heat capacity Cy and MCE. The 
vertical dashed line represents a cut at T = 0.08. Right Panel : the magnetic susceptibility on 
this cut. The inset displays the average magnetization revealing the m = 1/3 plateau. T and h 
are measured in units of Ji. 



4. The easy-axis J1-J2 Heisenberg model on a triangular lattice 

The metallic layered silver nickelate 2H-AgNi02 has been argued |13j to offer realisation of a 
spin-1 easy-axis AF on a triangular lattice with competing AF J2 interactions 

H = J,J2 S-S, + S.-S, -Dj2{Stf -hY^St (2) 

We have investigated the high field properties of Eq. [2] with ratio of parameters Ji = 1.0, 
J2 = 0.25 and D = 0.06 taken from preliminary fits to experiment [Hj. For these parameters, 
the classical ground state of Eq. [2] is a two-sublattice collinear Neel state with q = (0,7r). A 
natural expectation would be that, in applied magnetic field, this would undergo a l*** order 
spin flop transition into a canted two-sublattice state with the same wave number. However our 
results support a very different scenario. 

In Fig. [5]- (a) we show results of classical MC simulations of Eq. O compared with a mean 
field theory (MET) of the 2-sublattice spin-flop transition, and a 4-sublattice MET. The 2- 
sublattice theory predicts a first order spin fiop transition at h = 1.1. However a spin wave 




Figure 5. (Colour online) Left Panel : Ground state energies of Eq. [2]in 2-sublattice (dashed 
line) and 4-sublattice mean field approximations (blue line), together with results of classical 
MC simulation at T = 1 x 10~^ (red points). Simulation results follow the 4-sublattice solution 
for h > 0.12, avoiding the 2-sublattice spin-flop for h = 1.1. Right Panel : The quadrupole 
moment Q, showing a continuous transition into a state with transverse order. 



with q = (vr, ^^^) becomes soft at the much lower field of h = 0.12, precipitating a 2""^ 
order transition into a 4-sublattice "supersolid" state with both transverse order and broken 
translational symmetry. We can identify this transition using the quadrupole moment 

Q = ^Ql..y. + Qly Q.^-y^ = ^ E(^f )' - (Sff Q^y = ^Il '^sfsf (3) 

i i 

shown in Fig. [5]-(b). We have checked in interacting spin wave theory that this "supersolid" 
scenario also holds for quantum spins. Our preliminary results for the global h-T phase diagram 
also suggest that this supersolid state is superseded by cascade of phase transitions at higher 
field. These will be reported elsewhere. 
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